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The inherent asymmetry of the electric transport in graphene is attributed to Klein tunneling 
across barriers defined by pn-interfaces between positively and negatively charged regions. By 
combining conductance and shot noise experiments we determine the main characteristics of the 
tunneling barrier (height and slope) in a high-quality suspended sample with Au/Cr/Au contacts. 

We observe an asymmetric resistance Rodd = 100 — 70 LI across the Dirac point of the suspended 
graphene at carrier density |ug| = 0.3 —4-10 11 cm -2 , while the Fano factor displays a non-monotonic 
asymmetry in the range F a d a ~ 0.03 — 0.1. Our findings agree with analytical calculations based on 
the Dirac equation with a trapezoidal barrier. Comparison between the model and the data yields 
the barrier height for tunneling, an estimate of the thickness of the pn-interface d < 20 nm, and the 
contact region doping corresponding to a Fermi level offset of ~ —18 meV. The strength of pinning 
of the Fermi level under the metallic contact is characterized in terms of the contact capacitance 
C c = 19 x 10 -f> F/cm 2 . Additionally, we show that the gate voltage corresponding to the Dirac 
point is given by the work function difference between the backgate material and graphene. 


I. INTRODUCTION 

Klein tunneling is one of the most spectacular ef¬ 
fects of relativistic quantum field theory described by 
the Dirac equation. This tunneling phenomenon, present 
even in the regime of impenetrable barriers, leads to 
peculiar transport properties of graphene. Klein tun¬ 
neling is the backbone of transport due to evanescent 
modes causing the observed pseudodiffusive behavior of 
ballistic graphene samples ,2 . The bimodal distribution 
of transmission eigenvalues in ballistic graphene coin¬ 
cides with a diffusive conductor, which results in shot 
noise non-distinguishable from diffusive mesoscopic con¬ 
ductors. Furthermore, the evanescent modes lead to a 
minimum conductivity of in the ballistic regime . 
Evidence of these Klein tunneling phenomena have been 
obtained from observations of charge transport and shot 
noise in a graphene sheet with ballistic characteristics 1-1 ’. 

The most commonly employed assumption in the 
analysis of the conductance and shot noise of ballistic 
graphene has been to consider the carbon layer under¬ 
neath the electrodes as strongly doped, and this can be 
modelled using a rectangular electrostatic potential ’ ’ . 
In reality, this assumption suffers of severe limitations as 
in a real device the charge density varies continuously and 
the rate of change is governed by the screening length. 
Various theoretical models for finite-slope potentials have 
been analyzed for pn-interfaces in graphene “ . All the 
models have predicted asymmetry in transport prop¬ 
erties with respect to the gate voltage, i.e. whether 
the charge carriers are electrons or holes. In recent 
experiments, such asymmetry has been observed “ ’. 
In the ballistic regime, this asymmetry is attributed 
to the Klein tunneling “ while scattering by charged 
impurities ’ plays also a role in the diffusive regime. 
Furthermore, evidence of Klein tunneling has been re¬ 


ported in conductance experiments in confined geome¬ 
tries displaying phase-coherent and double-junction in¬ 
terference effects 18,19 . Sharp pn-interfaces have also been 
achieved in non-suspended samples using air-bridge type 
gates 20,2 . A full understanding of contact issues is of 
vital importance for the development of novel electri¬ 
cal components using graphene and other 2-dimensional 
materials . In particular, detailed understanding of pn- 
interfaces is critical for optoelectronics components '. 

The asymptotic carrier transport in Klein tunneling is 
bound to be affected by the strong influence of the metal 
contacts on graphene. A simple contact model was for¬ 
mulated by Giovannetti et air who also performed DFT 
calculations concerning the involved work functions. In 
this paper, we generalize this model to include the ef¬ 
fects of the applied backgate voltage, and we combine 
the resulting model with tunneling calculations based 
on the Dirac equation in order to obtain a comprehen¬ 
sive transport model for analyzing electrical conduction 
in a ballistic, suspended graphene sample. We employ 
a trapezoidal form for the tunneling barrier which we 
show how to treat analytically 2 1 . By using conductance 
and shot noise experiments performed on a high-quality 
suspended graphene sample, we can determine the bar¬ 
rier parameters and their relation to the doping of the 
graphene by a metallic contact. In comparison with DFT 
calculations , we find a semiquantitative agreement for 
the graphene-modified metal work functions as well as for 
the distance between the charge separation layers which 
govern the contact capacitance between the metal and 
graphene. 

The experiment and the agreement with the theoret¬ 
ical model confirm the existence of Klein tunneling in 
graphene. Our method works well even in the situa¬ 
tion in which the work functions of the contact metal 
and graphene differ by a relatively small amount (tens of 


2 


meV). As such, our results suggests a novel method to 
find the work function of materials. Note that the work 
function difference between two metals is not measurable 
directly: the standard way for its determination is the 
use of Kelvin probe force microscopy, where the electri¬ 
cal capacitance between the metal and a probe is varied 
in order to induce a measurable AC current. Our results 
demonstrate that there exists a gating effect in the posi¬ 
tion of the Dirac point of the suspended graphene due to 
the work function of the backgate. Thus, by using a ma¬ 
terial as backgate for graphene and measuring the gate 
voltage corresponding to the Dirac point one can get a 
simple DC measurement of the work function. 

The paper is organized as follows: after the present in¬ 
troduction as Sect. I, we present a complete theoretical 
treatment of the problem of gated suspended graphene 
with metallic contacts in Sect. II. The structure of our 
samples are detailed in Sect. Ill together with the em¬ 
ployed methods for shot noise measurements. Our exper¬ 
imental results and their analysis are presented in Sect. 
IV. The implications of the results are discussed in Sect. 
V jointly with a comparison to other works. 


II. THEORETICAL BACKGROUND 

A. Electrochemical model for suspended graphene 
samples 

When a graphene sheet is brought in contact with a 
metal, electrons will flow between them in order to equi¬ 
librate the Fermi level. This effect has been analyzed 
in detail in Refs. 24 and 26. We proceed beyond these 
works and introduce a complete electrochemical model 
for a suspended graphene sample with metallic contact 
electrodes. Our model takes into account consistently the 
effect of the backgate on the doping of the graphene un¬ 
der the metal. We show that this effect can be neglected 
only if the difference between the work function of the 
metallic contact electrode and the graphene layer is large 
enough (i.e. in the limit of high barrier, defined below 
Eq. (17)). If this is not the case, the contact-region dop¬ 
ing acquired from the backgate voltage has to be included 
in the transport calculations. On top of the electrostatic 
contributions, our model indicates that the work function 
difference between the backgate and graphene enters the 
equilibrium charge density, which, in particular, leads to 
a small shift in the Dirac point of the suspended part of 
the graphene. Using electrochemical equilibrium condi¬ 
tions, we derive analytical equations that incorporate all 
these effects. 

In Fig. 1 a) we present a schematic sideview of the 
sample, showing the region of contact with the metal 
and the graphene as well as the backgate. For clarity, 
we have exaggerated the size of some components, so the 
figure is not to scale. On the sample chip, the graphene 
sheet is placed at a distance of da = 300 nm from the 
backgate. The sheet is supported by a Si02 insulating 


layer ( e r = 3.9), which is partially etched away under 
the graphene. This results in a vacuum gap of height 
d V ac = 150 nm. Further details on the sample and the 
metallic contacts are found in Sect. III. 



FIG. 1. (Color online) Schematic of the sample and of the 
electrochemical potentials, a) Section perpendicular to the 
graphene sheet (pink), showing the metal M (grey) on top, 
a contact layer (yellow) of width d c , the support dielectric 
(beige), and the back gate (grey) g , at a distance da under 
the graphene layer. As seen in the cross section, we assume 
that the underetching of SiCb below the graphene can be re¬ 
garded as small, b) Schematic of the electrochemical poten¬ 
tials for the graphene under the metal (along the line M-g). c) 
Schematic of the electrochemical potential for the suspended 
graphene (along the line G-g). For definition of other symbols, 
see text. 

In our sample structure, the capacitance per unit area 
C g M between the graphene and the backgate in the sup¬ 
port region (graphene under the metal) is estimated from 
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the regular parallel plate formula 


C g M = 


eoe r 


for charge transfer between the metal and graphene is the 
difference between the work functions of the graphene un¬ 
der the metal and the work function of the metal, which 
is described by 


which yields C g m ~ 1.2 x 10 -8 F/cm 2 using the above 
mentioned values for do and e r . In the suspended re¬ 
gion, the graphene capacitance per unit area against the 
backgate C g c can be calculated using the formula for two 
capacitors in series: a vacuum capacitor with a plate sep¬ 
aration of d vac and a capacitor with the spacing d& — d vac 
filled with dielectric material having e r . This results in 
the capacitance per unit area 

C « G = , 77, _? w-i_ n - ( 2 ) 

ug T (uq d vac ) (e r 1) 


X = W G + A c - Wm, (4) 

where Wm is the work function of the metallic contact 
material and Wq + A c denotes the modified work func¬ 
tion of the graphene under the contact. Another electro¬ 
static potential U g m = en g M/C g M is established across 
the capacitance C g m, with a surface charge en g m on the 
gate. The total particle-number surface density «m hr 
the graphene layer in the contact region is therefore 

riM = n g M + n c . (5) 


In our calculations the actual value for the suspended 
part capacitance is taken from our Fabry-Perot measure¬ 
ments C g G = 4.7 x 10 -9 F/cm 215 , which agrees well with 
the above theoretical value. 

Next, we analyze the electrochemical potentials that 
appear in this experimental setting. We take two cuts 
through Fig. 1 a), one across the dashed line M-g, and 
the other across the dashed line G-g, and we represent 
the spatial variation of the electrochemical potential in 
Fig. 1 b) and Fig. 1 c), respectively. The work function of 
the pristine graphene is denoted by Wq. At the contact 
with the metal, this work function may be modified by a 
small shift, see Ref. 24. The work function of the contact 
metal on top of the graphene is denoted by Wm, while the 
work function of the backgate is W g . Note that accord¬ 
ing to Volta’s rule, the contact potential at the end of a 
circuit is determined only by the work functions of the 
circuit elements at the end; therefore no other work func¬ 
tion, for example corresponding to various other metals 
along the measurement chain, can enter in this problem. 
The difference between the work function of the gate and 
that of graphene is denoted by eVbirac = W g — Wq- this 
quantity will turn out to be the shift of the Dirac point 
of the suspended region of graphene due to the back- 
gate work function. Following the common practice in 
graphene research, gate voltages in the following formu¬ 
las will be measured with respect to the Dirac point, with 
the corresponding shifts defined as 

SVg = V g ~ V D irac- (3) 


The first equilibrium condition is obtained from the 
condition that the difference between the Fermi level of 
the metal and that of the gate equals eV g . This condi¬ 
tion does not formally involve the characteristic density 
of states of graphene. Hence, as shown by the circuit 
schematics below the Fermi level diagram in Fig. 1 b), 
it can be regarded as a pure electrostatic condition. It 
states 

e5V g = x + eU gM - eU c . (6) 

The second equation governing the equilibrium in¬ 
volves the intrinsic properties of graphene, and it can 
be obtained by using the condition that the Fermi levels 
of the metal and the graphene under the metal coincide: 

AI?fm + eU c = x- (7) 

In the graphene under the metal, where the linear 
graphene bands are supposed to persist, the relation be¬ 
tween the number of negatively-charged carriers per unit 
area um and the shift in the Fermi level A£fm is given 

by 


AE fm = hv F sgn [n M ]\A|n M |, (8) 

where vp = l.lx 10 6 m/s is the Fermi velocity. It is useful 
to introduce a constant £f relating the Fermi speed and 
the fundamental constants H and e: 


Cf = 



(9) 


We first solve the problem of finding the surface charge 
distribution for the electrochemical potentials presented 
in Fig. lb). At the contact between the metal and the 
graphene, the electrons will move from the electrode with 
the lower work function into the electrode with the higher 
work function. As a result, a surface charge distribu¬ 
tion en c will appear in the contact region, producing an 
electrostatic potential U c = en c /C c across the contact 
capacitance C c , the magnitude of which reflects the spa¬ 
tial variation of the surface charge. According to DFT 
calculations, C c ~ 10~ 5 F/cm 2 . A relevant parameter 


We propose to call this quantity Fermi electric flux. This 
constant is related to the concept of quantum capacitance 
(for graphene, see Ref. 27) and to the fine structure 
constant of graphene, as detailed in Appendix A. The 
Fermi electric flux determines the energy shifts produced 
by graphene as it is inserted in to an electrical circuit. For 
v F = 1.1 x 10 6 m/s, the equation yields Cf = 1.283 x 10~ 7 
V-cm. 

Combining now Eqs. (5-8) above, and choosing the 
proper physical solution, we obtain the final result for 
the shift of the energy level of graphene under the metal, 
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AE fm = sgn 


SV a 


xC c 

eC, 


sMj 



+ Cf Cc 


X + 


c, 


_gM 

C c 


eSV a 


( 10 ) 


An experimentally relevant limit for Eq. (10) is the case 
of a very large contact capacitance C c : if the material 
parameters are such that ( C c + Cm)Cf ^ Ixl and the 
charge induced by the gate voltage is relatively small, 
then eC g M\5V g \ <C C c |x|, and we obtain AEfm ~ X from 
Eq. (10). In this situation, the large contact capacitance 
locks the position of the Dirac point of the graphene un¬ 
der the metal to a value AEfm set by the parameter x- 

In our experimental setup, in fact, x remains relatively 
small, which makes the quantity eCM$V g /C c comparable 
with % at large values of SV g . Consequently, the Fermi 
level can even reach the Dirac point of the graphene un¬ 
der the metal AEfm = 0. As the latter condition does 
require substantial gate voltages, we could not reach this 
regime in our experiment. 

Next we turn to the suspended part of the graphene. 
In this region, the charge accumulated on the surface of 
the backgate n g G , which produces a voltage drop U g Q = 


zn g G/C g G, is exactly compensated by the charges on the 
graphene side having the particle-number surface density 
riG, 


n g = n g G- 

(11) 

Then, the condition that the electrochemical potential 
difference between the Fermi level of graphene and the 
Fermi level of the backgate is eV g reads 

eSVg = eU g G + AEfg, 

(12) 

where 


AEfg = hv F sgn [n G ]\A|nG|- 

(13) 


in parallel to Eq. (8). From Eqs. (11-13), we find for the 
Fermi level shift in the suspended part: 


AE fg = sgn [5V g \ 



eC gG (*\6V g 


(14) 


Note that AEfg = 0 corresponds to SV g = 0, i.e. 
V g = Voirac is indeed the Dirac point of the suspended 
graphene. For AEfg, typical experimental conditions 
fulfil e\SV g \ |x| ^ CgcCli which allows one to approx¬ 
imate AEfg ~ sgn [ 5V g \ Cf \J eC g c \5V g \. This limit was 
assumed in the previous analysis by Sonin . 

The difference of the Fermi energies of the suspended 
graphene and the graphene under the contact results in 
an electrostatic potential difference Vo between the free¬ 
standing region of graphene and the one under the metal: 


eVo = AEfm — AEfg- (15) 


This leads to a potential barrier for the electrons travel¬ 
ing along the ^-direction, from one metallic contact to the 
other, as detailed in the next subsection. Since transmis¬ 
sion through the potential barrier is foremost sensitive to 
the slope of the barrier at the charge neutrality point, we 
have adopted a trapezoidal barrier shape where the slope 
is constant on both sides of the flat top of the barrier. 
This trapezoidal form is tractable using analytical calcu¬ 
lations, and it is expected to catch the basic features of 
the tunneling transport problem. 


B. Klein tunneling: conductance and Fano factor 

The asymmetry in electrical transport across a basic, 
back-gated graphene device is related to many aspects 
of the sample, its quality, and biasing conditions. It de¬ 
pends on the nature of transport, whether it is ballistic 
or diffusive, on the presence of interfaces between p- and 
n-doped regions, on the steepness of the slopes of the elec¬ 
trostatic potential barrier, on the doping due to the con¬ 
tacts, and on the coherence of the transport between the 
reflecting interfaces (pn-type of interface or two unipolar 
regions with different doping, either p and p ’ or n and 
n\ yielding pp’- and rm’-interfaces). The type of junc¬ 
tion is determined by the sign of the Fermi level shifts: 
if AEfm > 0 and AEfg > 0 we have a nn' junction at 
the left contact, if AEfm < 0 and AEfg > 0 we have a 
pn junction, if AEfm > 0 and AEfg < 0 we have a np 
junction, and if AEfm < 0 and AEfg < 0 we have a pp' 
junction. The same notions apply to the right contact 
provided that the regions are considered from right to 
left, instead of the left-to-right direction used for the left 
contact above. 

Our experimental data deal with suspended graphene 
with a mean free path on the order of the sample length. 
Consequently, we will consider foremost Klein tunnel- 
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ing in the ballistic regime and neglect the influence of 
disorder* . Unlike the early work discussed above, we 
take into account the weak doping of contact regions 
when using Au/Cr/Au leads. We assume symmetric con¬ 
tacts although in reality there is always slight asymmetry. 
In our context, this assumption means that the values of 
the work function difference x are equal at both contacts. 
The asymmetry of the contact resistance can be neglected 
as the contact resistance is found to be insignificant in 
the analysis. The use of moderate voltages and cryogenic 
temperatures also guarantees negligible role of acoustic 
phonons . 

Below, we derive the formulae for the transmission and 
reflection coefficients, as well as the conductance and 
Fano factor, for the tunneling geometry given in Fig. 
2. The figure shows variation of the electrostatic po¬ 
tential (see Fig. 1) in the ^-direction along the whole 
graphene sheet including the end contacts. We consider 
only propagating modes; this approximation is justified 
by the criterion that propagation through the evanescent 
modes can be neglected if the bias voltages are clearly 
above hvpI&L , where L is the length of the sample. For 
a typical sample of length L = 1 /im, hvp /eE = 0.7 mV, 
well below the bias voltages used in our experiments for 
shot noise (10 — 70 mV). 


the contact area x < —d and L + d < x with zero (or 
small bias) potential, corresponding to the region M of 
graphene fully covered by the contact metal. Dirac elec¬ 
trons propagating along the ^-direction may thus experi¬ 
ence Klein tunneling under this barrier. In the left slope 
region — d < x < 0, the potential is written in the form 
V(x) = V 0 +sgn [Vo] a'^x. The parameter a sets the ab¬ 
solute value of the voltage slope. In general, the barrier 
can have either a positive or a negative slope, depending 
on the sign of Vo. For clarity, we will give the explicit an¬ 
alytical forms of the solutions of the Dirac equation only 
for the case of positive slope, corresponding to the left 
side of the barrier in Fig. 2. The solutions for negative 
slope can be obtained by the same procedure. 

All the energies are measured from the Dirac point of 
the graphene under the metal, with the standard con¬ 
vention that AEfm is positive if the Fermi level is above 
the Dirac point and negative otherwise, see Eqs. (8) and 
(10). A similar convention is used for AE F q, see Eq. 
(13-14). 

The massless Dirac equation for a graphene sheet can 
be written as 


. a .d 

ox oy 


^(x,y) = /C(a:)4'(a;,y), 


(16) 



FIG. 2. (Color online) Schematic view of a trapezoidal nn’n- 
barrier configuration across the sample. Vo denotes the differ¬ 
ence in the electrostatic potential between the contact regions 
and the center. The solid red line shows variation of the elec¬ 
trostatic potential along the graphene sheet (see also Fig. 1); 
the kinks of the red line separate the five different spatial re¬ 
gions in which the Dirac equation is solved. The slope of V ( x ) 
in the range —d < x < 0 is ahv-p/e. 


The use of the Dirac equation with a trapezoidal bar¬ 
rier is the simplest approximation for transport in a two- 
lead graphene sample with a pair of pn-interfaces due to 
spatially varying charge doping. The electrostatic trape¬ 
zoidal barrier (of structure nn’n) in Fig. 2 is comprised 
of five distinct regions: the middle part 0 < x < E, 
corresponding to the region G from the previous subsec¬ 
tion, where there is a step in electrostatic potential Vo 
(difference in the potential between the contact regions 
and the center), two sloped regions — d < x < 0 and 
L < x < L + d in the vicinity of the graphene/metal 
interfaces, where the electric field is finite, and finally 


where for particles near the Fermi level KL{x) = [AEfm — 
eV(x)\/Hvp, and the upper and lower components of the 
spinor 'F(a;, y) are denoted by ip+[x,y) and if>-(x,y ) re¬ 
spectively. Note that the electrons at the Fermi level will 
have energy AEfm hr any of the regions M, G, or in the 
slope region. We can solve the Dirac equation in Eq. 
(16) by the method of separation of variables, writing 
V’±( a b2/) = l/kt ( x ) sxp(ik y y), with the current along the 
x direction j x = ev F {\il)+\ 2 — l'0-| 2 ) normalized to ±er>F- 
The height of the energy barrier eVo = AEfm — AEfg 
defines a step eVo/Hvp in the momentum of the Dirac 
electrons as they cross the barrier. The influence of this 
change in momentum can be characterized compactly us¬ 
ing a dimensionless parameter, the impact parameter , de¬ 
fined as 


eVp 

hv-py/a 


(17) 


The impact parameter will play an essential role in our 
calculations below. Note that the earlier theoretical 
treatment by Sonin assumed a large impact parame¬ 
ter |po| 3> 1. As this approximation is not valid for our 
present experiment, a finite po needs to be taken into 
account in the theory. 

Assuming that the thickness of the sloped region d is 
independent of the gate voltage, allows us to express the 
absolute value of the slope of the barrier in terms of the 
V^-dependent quantities, AEfm and AEfGj which yields 


a = 


vi 

Cf d 


\Vo\- 


(18) 
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Likewise, we obtain for the impact parameter 



Note that the impact parameter increases with the ab¬ 
solute value of the slope, as evident from Eqs. (17) and 
(18). 


depending on weather the slope is positive or negative. 
Here x$ is the crossing point, defined as the position 
where the kernel is nullified, and it is given by Xq = 
sgn [Vo] A.Efg/clHvf. The solutions have been found by 
Sauter (see Ref. 25); for positive slope, with the substi¬ 
tution £(x) = \/a{x — xq) we have 


1. Top of the harrier, x > 0 (region G) 

In this region we consider an electron moving to the 
right, j x = evp. The Dirac equation can be solved, with 


M x )\-d<x<o = C^iax^ky/y/a) + C 2 G* (£(&), k v /y/a) 

lM*)|— rf < x <o = C 1 G(ax),k y /V^)+C 2 F*(^),k v /V^). 


K(x\x > 0) = — AE fg , (20) 

hv f 

for both positive and negative AEfg, 


Note also that the dimension of a is [length] -2 , while 
that of K(x) and k y is [length] -1 ; hence the quantity 
ky/y/a, that enters the hypergeometric functions, is adi- 
mensional. 


( VhO) 

V 


^ x>0 


1 fex,FG=l -iky 

2 ' 2 k FG 

_ 1 _i_ ^x,FG 

2 2fcpG 



±ik XfFG x 


( 21 ) 

where the ± sign denotes sgn[AE F G]- The absolute value 
of the wavevector is 


fcFG 


|AF F g| 
hv F 


( 22 ) 


The functions F and G are defined through the Rum¬ 
mer confluent hypergeometric function M = \F\, 


F(£,k) = e~^ 2/2 M 



(24) 


and k x ,FG = y^FQ — ky. Note that k y remains un¬ 
changed over all region crossings. 


2. Slope of the barrier, —d < x < 0 


and 


G(£, k) = —n£,e 


-i£ 


'I 2 m( 


1 - 


4 2 




(25) 


In this region we have 

K,(x | — d < x < 0) = —sgn [Vq] a(x — xq), (23) 


From the boundary conditions at x = 0 we obtain 


C i 

g 2 


kFG + k x ,FG + *sgn[AF F G]fc y 


2\A'FGk; 


F*(?(0),fc y /v^) 


&FG — k x ,FG — *sgn[AF F G]fcy n * 


i,FG 


2-\A'FGfcai,FG 


G*im, ky/y/a), 


-kFG + k x ,FG+^gn[AE F G}ky F{mky/Va) _ fcFG + k x , F G + »Sgn [AE F G]fcy 
2-\AfG& 2;,FG 2yfc F Gfca:,FG 


(26) 

(27) 


Importantly, a consequence of the normalization \i/>+\ 2 — 
|^>_| 2 = 1 is the fact that |Ci| 2 — IC 2 I 2 = 1 and 
|F(£, k)| 2 — |G(£, k )| 2 = 1, which can be verified explicitly 
using the explicit expressions above. 

3. Graphene under the metal x < —d (region M) 

In the calculation of Sonin i5 , this region with x < —d 
was disregarded on the basis of the assumption p\ 


1). It turns out that this condition is not valid in our 
experiment (cf. Fig. 5), and the behavior in the region 
x < —d has to be taken into account. In this region, 

IC(x | — d < x < 0) = -—AFfm- (28) 

hv f 

The absolute value of the total momentum at the Fermi 
level is then 

^FM = T-|AFfm|- (29) 

hv p 
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The corresponding momentum in the x-direction equals 


to k.. 


.FM 


= Jk 


‘FM 


— ky , and the wave is a superposition 


of a reflected and a transmitted component, 


( V>+0) 

V 


x<—d 


1 

t 


1 fca,FM=l -^y 

2 2fepM 

_ 1 fcl.FM ±ik y 
2 2/cfm 



± ik x 


Here r and t are the reflection and transmission ampli¬ 
tudes, ± is the sign of sgn[AFpM]) and =F in the last 
exponent is —sgnjAFpM]- 

Note also that in Eq. (30) the reflected component 
is normalized to the current in the x-direction, equal to 
—euF, while the transmitted component is normalized to 
+euF- However, one can explicitly check that the overall 
normalization of U/(a;,y) is to +ev F , provided that |?’| 2 + 
|f| 2 = 1. This ensures that the normalization is the same 


( 1 — fex,FM ±iky 

2 2/cfm 

_ 1 I —k Xy FM±ik y 

2 2 k FM 



■F^ 3 :,FM® 


(30) 


for all the three regions. To understand intuitively how 
this is realized, note that the admixture of the reflection 
component in H/(x,y) of Eq. (30) is compensated by an 
increase in the component propagating to the right, since 
\t\ becomes subunitary. 

By imposing the condition of continuity of the wave 
function at x = —d, after some algebra we obtain the 
complex transmission amplitude 


1 FMrf _ _~ k FM - k x , FM + isgn[AE FM ]ky 


-e 


2\/^FMfci,FM 


[CiF(£(-d), ky/y/a) + C 2 G*(tt-d), ky/yTa)} 


feFM - fc*,FM + *Sgn AFfm fcj,. . . , / n „ „*/*, « , , n , 

- = - u -{CiG(t,{-d) 1 ky/yJa) + C 2 F (£(-d),ky/y/a)]. 

2-*/ kfm%,fm 


(31) 


C. Conductance and Fano factor for the whole 
barrier 


To calculate the total transmission through the bar¬ 
rier, we employ incoherent addition of the transmission 
coefficients. Usually phase coherence is more sensitive to 
disorder than reflection and transmission, and we address 
the case when the former is destroyed but the latter is not 
affected by disorder (ballistic regime). This assumption 
is supported by the experimental fact that the Fabry- 
Perot resonances are found to be weak. Furthermore, 
when destroying the Fabry-Perot resonances fully by an 
applied bias, the overall conductance does not change 
much. Hence, we consider the incoherent treatment of 
transmission probabilities well justified in our analysis. 

In general, for the case of incoherent tunneling through 
a symmetric barrier with equal transmission T for the 
left and right slopes, the total probability of transmission 
through the barrier is (see Appendix B) 


Ttot 


1 

2T- 1 - 1' 


(32) 


By using the Landauer-Biittiker formalism, we obtain 
the conductance and the Fano factor as sums over the 
transmission coefficients and their quadratic values. Each 
quantized value of k y corresponds to a conduction chan¬ 
nel, over which the summation of transmission coeffi¬ 
cients has to be performed in order to obtain the total 


conductance (shot noise) from the conductance per chan¬ 
nel (shot noise per channel). Thus, for a sample with a 
given level of contact doping A£fm, one can calculate 
the conductance a and the Fano factor F as a function 
of gate voltage. In the limit of large number of channels, 
a and F can be written as 


4e 2 IU 

|.min{fepM,fepG} 


7 xh 

L 

dky'l tot 

1 

y'min{fcpM,fcFGl 


F = 1 - 

/ 

dk y T tot . 


(33) 

(34) 


The conductance and the Fano factor will depend on 
V g through the dependence of the quantities /cfGi ^fm, 
A£fg, A£fm obtained previously. If the contact resis¬ 
tance between the graphene and the metal is neglected, 
then the input to Eq. (32) is given by T = \t] 2 , where 
|f| 2 = \t(k y )\ 2 is the transmission coefficient calculated 
using Eq. (31). 

If a finite contact resistance exists, then a finite trans¬ 
mission probability 0 < T c < 1 should be included in the 
value of T in Eq. (32). By applying again the result of 
Appendix B, we may write 


T = 


Tc\t \ 2 

T c + \t\ 2 — T c \t\ 2 ' 


(35) 


The inclusion of contact resistance has a strong influence 
on the shot noise, and the calculated Fano factor becomes 























quickly larger than the measured value when the contact 
transmission is lowered from one. 

The limits of integration in Eqs. (33) and (34) are set 
by the condition that the wave vector k x is a real number, 
in other words the electron is not in a bound state but 
propagates to infinity. The condition k y < kpG comes 
from the top of the barrier region, while the condition 
k y < kpM comes from the region x < —d. 

Interestingly, even though our potential profile with 
five distinct regions as shown in Fig. 2, the entire model 
has only two essential fitting parameters, namely po = 
eVo/hvpy/a = sgn [Vo] yfad and the barrier slope a. The 
impact parameter contains the information on the dop¬ 
ing via Eqs. (8), (14), and (15) and it enters in the 
upper limit of the integrals in Eqs. (33-34), since fcpM = 
Iv^Po + AEfg/Hvf\- In fact, the Fano factor in our cal¬ 
culation is fully determined by the value of po> while the 
conductance integrals need also the value of the barrier 
slope for their evaluation. 


III. SAMPLE AND SHOT NOISE METHODS 

The measured suspended sample was manufactured 
using standard PMMA-based e-beam lithographic tech¬ 
niques on a graphene piece exfoliated onto Si/Si 02 chip 
(see Fig. 3). The dashed white line with two circles in 
the figure indicates the positions at which one obtains 
the cross sectional view depicted in Fig. la; the cir¬ 
cles correspond to the spots at which the electrochemi¬ 
cal potential profiles (along vertical direction) are drawn 
in Figs, lb and lc. The graphene sample was con¬ 
tacted with Au/Cr/Au leads of thickness 5/7/50 nm, 
after which roughly 135 nm of the sacrificial silicon diox¬ 
ide was etched away using hydrofluoric acid (HF) fol¬ 
lowing the methods discussed in Refs. 30 and 31. Ra¬ 
man spectroscopy was employed to verify the single layer 
structure of the graphene sample. Current annealing at 
V g = 0 was used to enhance the mobility of the sample. 
We employed voltage bias around 1 V and a current of 
0.3 — 0.7 pA/pm in our cleaning process. The aspect ra¬ 
tio of the sample as determined before the experiments 
was W/L = 4.5 pm/l.\ pm ss 4.1. The capacitance 
CgG =47 aF/ pm 2 was determined using Fabry-Perot in¬ 
terference fringe measurements . The mobility of the 
sample pa was calculated using the charge carrier den¬ 
sity |nc| = \C g G,5Vg\/e (obtained from Eq. (14) in the 
limit of large 5V g ) and a{6V g ) = R~ 1 (5V g )L/W in the 
formula po = (er — &™in S )/ n G e i where is the mea¬ 

sured minimum conductivity. We find pc > 10 5 cm 2 /Vs 
near the Dirac point at n < 2.5 • 10 10 cm -2 . For the 
Fermi velocity we used the value vf = 1-1 • 10 6 m/s; note 
that owing to interaction effects at small charge density, 
the Fermi velocity can grow up to vf — 3 • 10 b m/s in our 
sample near the Dirac point ’. The assumption of sym¬ 
metric contact capacitances (within ±15 %) was verified 
from the inclination of the Fabry-Perot pattern (see Ref. 
15). 



FIG. 3. (Color online) False color scanning electron micro¬ 
graph taken of a suspended sample similar to the measured 
one (topview on y-x plane; x-axis is oriented vertically): yel¬ 
lowish areas denote Au/Cr/Au contacts while the greenish 
area in between depicts the suspended part of the graphene 
sheet. Wrinkles due to strain are visible in the sheet at room 
temperature. Charge is induced into the central graphene re¬ 
gion by the backgate (formed by the doped silicon substrate, 
not visible) and, consequently, a pair of junctions (abrupt 
changes in the charge density) are formed very close to the 
contacts. The upper white circle denotes a position where the 
schematic profiles of Fig. lb apply, while the profile at the 
lower circle is governed by Fig. lc. The scale is given by the 
white bar on the right. 


Differential conductance of the sample was measured 
using standard low-frequency lock-in techniques (around 
35 Hz); the same excitation was also employed in our 
differential shot noise measurements. The noise-signal 
from the sample was led via a circulator to a cryogenic 
low noise pre-amplifier having a bandwidth of 600-900 
MHz . The amplifier provided a gain of 15 dB and the 
signal was further amplified at room temperature by 80 
dB and band pass filtered before the Schottky diode de¬ 
tection. Small bonding pads of size 90 pm x 90 pm were 
employed in order to keep the shunting capacitance ~ 0.1 
pF negligible in our microwave noise measurements 5 ’. 
For calibration purposes, a microwave switch was used 
to select a tunnel junction as the noise source instead of 
the sample. For details of the calibration procedure, see 
Ref. 33. In our experiments, we employ the excess Fano 
factor obtained from the noise: F = (Si(Vi,)— 5/(0))/2e/, 
where we use the difference of the current noise spectrum 
Si between bias voltage V/ = 22 mV and V), = 0. The 
experiments were performed around 0.5 K in a Bluefors 
BF-LD250 dilution refrigerator. 


IV. MEASUREMENTS AND RESULTS 

Before presenting our data and their analysis, let us 
first summarize the free parameters that enter our theo¬ 
retical model. The capacitances are fixed by the geome¬ 
try and the charge dipole layer thickness Zd , which equals 
to Zd = 0.9 A according to the calculation of Ref. 24 for 
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graphene/gold interface; assuming Zd corresponds to d c 
in Fig. 1, this yields C c = 9.8 x 10 -6 F/cm 2 using the 
vacuum permittivity. There are three parameters that 
are fitted to the data: % = Wq + A c — Wm , AEfm, and 
the thickness d of the charge gradient interface. The over¬ 
all magnitude of the measured conductance sets % = —18 
meV in our analysis. Furthermore, the noise analysis in¬ 
dicates that d ~ 20 nm (see below). Hence, we are left 
with one (5F/-dependent fit parameter, AEfm, which de¬ 
scribes the electrostatic behavior of the contact regime 
with the gate voltage. 


The zero-bias conductance a as a function of V g is dis¬ 
played in Fig. 4. The Dirac point resides at Foirac = 
—0.2 V indicating negative dopants on the sample. How¬ 
ever, the asymmetry of the conductance and the shot 
noise suggest positive doping, which would correspond 
to Foirac > 0 V, in contrast to the above. In order to 
account for the ’’wrong” sign of the Dirac point location, 
we argue that the work function of the back gate has to 
play a role. The work function of doped silicon depends 
on the sign and amount of the dopants. At large nega¬ 
tive doping, W g — 4.4 eV for our background material 
with a negative dopant concentration of ~ 10 16 cm -2 . 
Since Wq = 4.6 eV according to Ref. 35, we obtain 
Foirac = W g — Wg = —0.2 V, which leads to SV g = 0 at 
V g = —0.2 V. Hence our data are consistent with doping 
induced by the work function of the backgate material. 


As in earlier works, we characterize our results in terms 
of the odd part of the resistance across the Dirac point, 
-Rodd = (R(+5Vg) — R(— SV g ))/2 with SV g counted from 
the Dirac point, we find R 0 dd = 100 — 70 at carrier 
density |jig| = 0.3 — 4 • 10 11 cm -2 , which is consistent 
with other experiments 1 “ . The conductance at large 
positive SVg is substantially smaller than at equal carrier 
density at SV g < 0 as expected in the presence of pn- 
interfaces. The overlaid curves in Fig. 4 represent the 
fits of the model for each gate voltage value using the 
Fermi level position AE FM (V g ) as the fitting parameter, 
with set values for the thickness d = 20 nm and for the 
work function difference \ = —18 meV; further input 
parameters of the model are specified in Sect. II A. Since 
the fitting is done separately at each point, there is no 
difference between the overlaid curve and the data. 



-15 -10 -5 0 5 10 15 

v v > 

FIG. 4. (Color online) Zero-bias conductance R~ x as a func¬ 
tion of the gate voltage V g ; V g > Foirac corresponds to a pnp 
structure, while V g > Fbirac represents data in a pp'p con¬ 
figuration. The overlaid trace depicts the curve which was 
employed to determine the Fermi level shift AEfm(F 9 ); this 
Fermi level shift is displayed vs. SV g = V g — Foirac = V g + 0.2 
V in Fig. 5. 

Fig. 5a depicts the Fermi level shift AEfm in the con¬ 
tact regions, obtained from the point-wise fitting of our 
model to the conductance data in Fig. 4. The Fermi level 
does not move much with SV g , which is due to a large 
contact capacitance C c when compared with the other 
capacitances in the system (see Eq. (10)). Neglecting 
the region around the Dirac point, AEfm changes almost 
linearly with gate voltage due to the additional charge in¬ 
duced by V g to the contact region. This is expected from 
our model, and by fitting Eq. (10) to our data in Fig. 5a 
we obtain % = —18 meV, C c = 1.2 x 10~ 5 F/cm 2 , and 
C g M = 6x 10 -9 F/cm 2 , close to the geometrical estimates 
discussed earlier. Taking into account the uncertainties 
in the involved capacitances, the asymptotic agreement 
between the experiment and the theoretical model can 
be considered as good. 

Additional effects are observed near the Dirac point 
<5F/ = 0. The doping AEfm becomes close to zero, and 
clearly our electrostatic model does not capture the be¬ 
havior completely. One shortcoming, for example, is that 
we do not include the renormalization of Fermi velocity 
by interactions near the Dirac point, which would partly 
improve the agreement between Eq. (10) and the mea¬ 
sured data. Also near the Dirac point evanescent waves 
become important, but they can be neglected when work¬ 
ing at finite Ft, 2 Yivf l&L. Finally, another explanation 
could be the formation of charge puddles, which exist in 
suspended graphene even though their strength is sup¬ 
pressed compared with non-suspended samples”. These 
phenomena are beyond our present analysis. 

Fig. 5a displays our experimental results for AEfg 
and the difference A£ FM — AE F g ■ Once AE F m — AE F g 
is known, we may evaluate the impact parameter po = 
eVo/hvFy/a, which is depicted in Fig. 5b. We note that 
the value of the impact parameter is proportional to \fd 
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(when Vo is fixed, a oc 1/d), apart from small correc¬ 
tions. At large interface thickness, our data approaches 
the regime of validity of the analysis by Sonin . Accord¬ 
ing to Eq. (18), the 6V g dependence of the slope of a in 
our barrier is given by AEfm — AEfg when the interface 
thickness d is fixed. 
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formula in Eq. (34) for interface thicknesses d = 1, 20, 
and 100 nm between the doping sections dominated by 
the metal and the backgate. The dependence of the cal¬ 
culated Fano factor with the interface thickness is found 
to be rather weak below d = 20 nm. 

As with resistance, we have analyzed the odd part of 
the Fano factor, F odd = F(+SV g ) — F(—6V g ) , which 
would be independent of spurious, common contributions 
due to scattering on both sides of the Dirac point. For 
example at 5V g = 72 mV, we find an asymmetric maxi¬ 
mum at |fc| ~ 3 • 10 5 cm -1 above which F odd decreases 
smoothly from ~ 0.1 down to 0.03 at charge density of 
|no| = 4-10 11 cm -2 . The experimentally obtained value 
F odd ~ 0.1 and the slow decrease with \k\ correspond well 
to the theoretical behavior. 



FIG. 6. (Color online) Measured Fano factor (circles) vs gate 
voltage V g , and results from the contact model using parame¬ 
ters obtained from the conductance; the curves are calculated 
for thickness d = 1, 20, and 100 nm (blue, red, green, respec¬ 
tively) . 


FIG. 5. (Color online) (Upper frame) Position of the Fermi 
level AFfm in the contact region obtained from the conduc¬ 
tance analysis of Fig. 4 (orange symbols), the Fermi level 
position in the suspended part AEfg (blue symbols), and 
the difference AEfm — AEfg (red symbols). The solid line 
is the dependence obtained from Eq. (10) using x = ~ 18 
meV; the values of the contact capacitances are discussed in 
the text. (Lower frame) Variation of the impact parameter 
po with gate voltage obtained from the AAfm — AAfg data 
presented above; po < 0 corresponds to the pnp configuration. 
The behavior depends on the thickness d, the values of which 
are indicated in the figure. 

Last, we present our data on the Fano factor. We em¬ 
ployed a bias voltage 14 > 20 mV in our measurements 
in order to be well above the cross-over voltage between 
thermal and shot noise, as well as to reach clearly the 
regime where the incoherent summation of barrier trans¬ 
missions becomes valid (see Subsection IIC). The exper¬ 
imental data measured at Vj, = 22 mV are presented in 
Fig. 6 together with the curves given by our theoretical 


Our noise data clearly display asymmetry that is ex¬ 
pected for Klein tunneling. Our model yields Fano fac¬ 
tors that are in good agreement with the measured re¬ 
sults at 5V g > 0 (the case with pn-interfaces) when the 
interface thickness is set to 1 nm. Thick interfaces on the 
order of 100 nm do not fit the results whereas intermedi¬ 
ate thicknesses d ~ 20 nm and below fit reasonably well. 
Consequently, we conclude that d < 20 nm on the basis 
of our shot noise results. 

At 5V g < 0, our model predicts a clear drop in the Fano 
factor with decreasing SV g and then a subsequent grad¬ 
ual recovery with lowering electrochemical potential. The 
data displays similar features but the SV g scale appears 
much smaller. The calculated Fano factor at SV g < 0 
appears to be in agreement with the data on average . 
The data display a clear dip on the negative gate voltage 
side. Asymptotically, there is a clear difference in the cal¬ 
culated Fano factors, as in the experimental results. The 
agreement could be improved by adding a fitting param¬ 
eter to account for a small reduction of contact trans¬ 
parency, but here we have set the contact transparency 
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T = 1 in order to avoid an additional fit parameter. We 
believe that these differences are mostly caused by non 
linear screening effects which are known to be different at 
the Dirac point and away from it ’ . Such charge den¬ 

sity dependent screening is not included in our model. 


V. DISCUSSION 

In the experiments on graphene, it is common that the 
Fabry-Perot type interferences remain weak even though 
the sample is more or less ballistic 1 . This is accordance 
with the starting point of our analysis in which the in¬ 
terference effects are neglected. Indeed, phase coherence 
is satisfied over small distances on the order of 20 nm 
(i.e. over the thickness of the barrier) as assumed in our 
calculations. However, the coherence is almost lost over 
the distance between the barriers and interferences be¬ 
come negligible. Besides, our noise experiments were per¬ 
formed at finite bias, above the regime where Fabry-Perot 
oscillations were observed, which may strengthen the 
tendency towards incoherent behavior; in fact, electron- 
electron interaction effects in graphene have been found 
to be 100 larger in graphene than in regular metallic 
systems 1 . 

The work function difference including the chemical 
interaction A c was found to be related to the over¬ 
all magnitude of the sample conductance, this did set 
X = Wg — Wm + A c = —18 meV. In our notation, x < 0 
corresponds to positive doping if AUfm is solely governed 
by x- For gold contacts (Wm = 4.7 eV), DFT calcula¬ 
tions predict x — —0.1... — 0.2 eV 24 ’ 39 , but one has to 
keep in mind that these calculations yield the work func¬ 
tion of pure gold with an error on the order of 0.2 eV. 
Experiments suggest that x = —0.35 eV for pure gold 1 . 
Furthermore, Cr has a work function of Wm = 4.5 eV and 
Cr/Au contacts have been shown to yield a work func¬ 
tion of 4.3 eV for graphene under the contact 11 , which 
leads to an estimate x — —0.2 eV for Cr/Au contacts. 
The results of Ref. 42 suggest x = —0-1 eV for Au/Cr 
contact, although the authors discuss that the contact 
doping might be close to zero, which would be in agree¬ 
ment with our results. However, the exact comparison 
with Ref. 42 is problematic because of the differences 
in the contact structure: we evaporate first 5 nm of Au 
before laying down 7 nm of Cr. As discussed below, the 
vacancy creation by the first evaporated metal layer is of 
major concern. 

The actual electrical contacting is further complicated 
by the reactivity of metal atoms on top of graphene 1 ’’ 1 . 
Using in situ X-ray photoelectron spectroscopy it has 
been demonstrated that sometimes the side-contacting 
picture may be misleading with real contacts 1 '. In our 
analysis, for simplicity, we need to assume a uniform 
graphene layer under the metallic contact, although it 
is possible that metals like Cr promote vacancy forma¬ 
tion and lead to creation of defects under the evaporated 
metal; similar defect formation has been reported due to 


deposition of gold atoms as well as in annealing stud¬ 
ies of metallic contacts 1 . Furthermore, charge transfer 
at the interface depends on the amount of oxygen and 
nitrogen on top of graphene as has been found in func¬ 
tionalization studies by Foley et al . 

Despite of possible defect formation, gold contacts are 
known to preserve the graphene cone structure under 
the contact • . Quantum capacitance measurements of 

graphene under the contact were performed in Ref. 50 to 
characterize the cone structures. The accuracy of these 
quantum capacitance determinations compared with the¬ 
ory is within a factor of two. These findings suggest a 
modification of the Fermi velocity under the contacts. An 
inclusion of such a modification would form a promising 
extension of our theory, but this was left for future. 

The parameter AEfm together with gate capacitance 
(yielding A7?fg) determines fully the product da which 
varies strongly with the gate voltage. Hence, our mea¬ 
surement imposes a constraint on the product da as a 
function of 5V g . When the thickness of the interface is 
fixed, the slope of the interface approaches quickly zero 
as AEfm — AAfg —l► 0. Altogether, the range of varia¬ 
tion of AEfm — AEfg is rather limited, which indicates 
weak doping by contacting metal as well as by the gate 
(CgQ is small for suspended devices). 

The same analysis as presented here can be performed 
using a constant slope a rather than a fixed d. We did 
such an analysis and, somewhat surprisingly, the numeri¬ 
cal simulations yielded similar predictions for the conduc¬ 
tivity and Fano factor, in particular for i? 0 dd and F 0( jd- 
This increases the confidence in the analysis of our re¬ 
sults, demonstrating that the extraction of relevant pa¬ 
rameters does not depend strongly on model-specific as¬ 
sumptions about the trapezoidal barrier in the regime of 
our data (i.e. at small interface thickness d). 

Screening influences the speed at which charge den¬ 
sity varies across borders between differently doped re¬ 
gions. Due to its peculiar inherent properties, screening 
in graphene may be strongly non-linear. According to 
theory " , the screening influences mostly the asymptotic 
behavior of the barrier, whether it is aU 1 / 2 or a; -1 , but 
not much the length scale of the rapid initial relaxation. 
Since we use trapezoidal shape and neglect the relax¬ 
ation behavior in the barrier altogether, we have chosen 
to work with fixed thickness, which is proportional to 
the average inverse slope of the charge relaxation at the 
pn-interface. One should keep in mind that our analysis 
is not reliable close to the Dirac point since there the 
slope itself weakly affects the conductance. Therefore, 
our fitting is mostly sensitive to the interfacial thickness 
far away from the Dirac point where the screening length 
becomes short. 

According to Ref. 37, the approximate thickness of the 
pn -interface could be in the range of 5 nm (see also Ref. 
51), which is consistent with our results d < 20 nm. Full 
numerical simulations on the pn-interface structure in a 
double-gated graphene structure have been performed in 
Ref. 52. Within the trapezoidal approximation, the cal- 
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culated slope of the potential profile in Ref. 52 yields an 
interface thickness of 30-40 nm. In our case, the leads 
may act as gates and, consequently, the pn-interfaces 
near the contacts remain sharp, in agreement with theo¬ 
retical estimates. Our interfacial width of d ~ 20 nm, on 
the other hand, is much smaller than found using scan¬ 
ning photocurrent microscopy on non-suspended samples 
fabricated on silicon dioxide ’ . 

Our primary fit parameter AI?fm takes into account 
the standard electrostatics in the contact region. Our 
analysis indicates a clear success of electrostatic analysis, 
and the results verify the role of large contact capacitance 
arising due to charge transfer between the contact metal 
and graphene. This leads to pinning of Fermi level at 
the contact, which is consistent with findings in Refs. 
53 and 54. Near the Dirac point, we find modifications 
from the standard electrostatic doping picture, which are 
presumably related with the neglect of proper screening 
treatment and nonuniformities in the charge distribution 
near the Dirac point. 

Finally, our analysis is based on a rather idealized theo¬ 
retical model. Additional effects can be included in to our 
theoretical model, for example the broadening of the den¬ 
sity of states in the graphene due to inhomogeneities and 
due to coupling to metal. Clearly, the inclusion of these 
effects would broaden the conductance and the Fano fac¬ 
tor characteristics, resulting in a better fit with our ex¬ 
perimental data. The disadvantage of this phenomeno¬ 
logical approach is that supplementary knowledge on ad¬ 
ditional parameters would be needed. We did check, 
however, what happens if a finite contact resistance is 
included using Eq. (35). From the simulations, we find 
that a significant discrepancy in the overall value of the 
Fano factor starts to appear for t c < 0.9: the Fano factor 
increases quickly over the entire gate voltage range and 
the fitting becomes difficult, no matter how large doping 
is used. This constrains the value of t c to ~ 0.95 — 1, 
but the improvements in the fitting achieved within this 
range are negligibly small. 

In conclusion, we have developed a practical trans¬ 
port model for analyzing transport in ballistic, suspended 
graphene samples. Our analysis shows how the doping 
in graphene depends on the difference between the work 
function of the metal and graphene as well as on the ap¬ 
plied gate voltage measured from the Dirac point which is 
given by the work function difference between the back- 
gate material and graphene. By combining conductance 
and shot noise experiments performed on a high-quality 
suspended graphene sample, we have determined all the 
relevant parameters which are involved in the electrostat¬ 
ics of the contact and in the Klein tunneling of graphene. 
When comparing with DFT calculations 1 , we find a 
semiquantitative agreement for the graphene-modified 
metal work functions as well as for the distance between 
the charge separation layers which govern the contact 
capacitance between the metal and graphene. The small 
charge layer separation (~ 1 A) leads to a large contact 
capacitance which is responsible for a rather weak tun- 


ability of the Fermi level position under the contact. 


ACKNOWLEDGEMENTS 

We acknowledge fruitful discussions with D. Cox, 
V. Falko, T. Heikkila, M. Katsnelson, M.-H. Liu, 
and A. De Sanctis. Our work was supported 
by the Academy of Finland (contract 250280, LTQ 
CoE) and by the Graphene Flagship project. This 
research project made use of the Aalto Univer¬ 
sity Cryohall infrastructure. S.R. and M.F.C ac¬ 
knowledge financial support from EPSRC (Grant 
no. EP/J000396/1, EP/K017160/1, EP/K010050/1, 
EPG036101/1, EP/M001024/1, EPM002438/1) and 
from Royal Society international Exchanges Scheme 
2012/R3 and 2013/R2. 


Appendix A: The Fermi electric flux 


The Fermi electric flux 


Cf = 



(Al) 


is related to the concept of quantum capacitance per unit 
area of graphene, 


< A2 > 


where n is the number of excess (negatively charged) car¬ 
riers per unit area, V(AE F ) is the density of states of 
graphene, and A E F is the energy of the Fermi level mea¬ 
sured from the Dirac point, 


A E F = hv FSgn [n] \/n\n\ = e^FSgn[n]\/H- (A3) 

Another useful formula is |n| = (AEf/cCf) 2 - With these 
notations one immediately obtains 

AE F = l -C q Cl, (A4) 


resembling the formula for the energy ^CV 2 of a capaci¬ 
tor C charged by a fixed voltage V. Therefore the Fermi 
electric flux is the electric flux through the plates of a ca¬ 
pacitor of unit area and unit distance between the plates, 
with capacitance equal to the quantum capacitance and 
charged to an electrostatic energy A E F . 

The Fermi electric flux is a useful constant also for 
studying the inductive properties of graphene. The dual¬ 
ity between the electric and magnetic properties is man¬ 
ifested in this case through the simple relation = 
where $0 = h/2e is the flux quantum. Then, 
the graphene kinetic inductance can be defined by regard¬ 
ing the Fermi energy e F as a magnetic energy associated 
with the magnetic flux S, F /v F , yielding 


^G.kin 


JL 

v^e F 


(A5) 
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FIG. 7. (Color online) Generic schematic for tunneling 
through two interfaces, separating a middle region ( b ) from 
the left ( l ) and the right (r) regions. 


in agreement with the known result for the kinetic induc¬ 
tance of graphene . 

Another connection can be made with the fine struc¬ 
ture constant, a = e 2 /Ane^hc ~ 1/137, which has an 
essential role in determining the optical properties of 
graphene A We obtain 


Cf 


evp 

A^/tteq ac 


(A6) 


Similarly, one can introduce a graphene fine structure 
constant qq = e 2 /Ane^hvF = ac/v f ~ 2, obtaining = 

e/Ay/Tre 0 a G - 


Appendix B: Conductance for the whole barrier 

Equation (32) is well known in mesoscopic physics and 
can be found in some textbooks, e.g. by S. Datta and 
by Yu. Nazarov & Ya. Blanter . It is derived by sum¬ 
ming over all the Fabry-Perot reflection and transmission 
processes and neglecting the interference term in the fi¬ 
nal result . Here we give an alternative, simpler proof, 
starting from the beginning with the assumption that the 
currents simply add up without any interference term. 

We consider the generic problem of incoherent tunnel¬ 
ing through two interfaces, interface 1 separating the left 
region from a middle region and interface 2 separating the 
middle region from the right region, see Fig. 7. Typically 
in this type of problems, the middle region is associated 
with a barrier, so we will denote it as such. The barrier 
has transmission coefficients T-\ on the left side and T 2 
on the right side; the corresponding reflection coefficients 
are denoted by i?i = 1 — T\ and R 2 = 1 — T 2 . 

We assume that a wave producing a current prop¬ 
agates from the left ( l ) region towards the barrier. Part 
of this wave will get reflected back to the left into j^), 
and part of it gets transmitted in the region under the 
barrier (6), with the current j^}. At the right side of the 
barrier, the incoming current //A gets transmitted to the 
region on the right of the barrier (r) as j^} and part of it 
is reflected back under the barrier as . The wave j^} 
travels backwards (to the left) towards the left side of the 
barrier, where it is partly reflected and partly transmit¬ 
ted, contributing to and j^} respectively. By using 
the fact that in the absence of interference the currents 


will just sum up, we get 


3$ = T 2 jW, 

(Bl) 


(B2) 

jib) = Tij W + Rij w t 

(B3) 

= Rij® + TxjW. 

(B4) 


Adding the first pair of equations Eqs. (B1-B2) yields the 
conservation law for the current at the (6) — (r) interface, 
j= jfW — jA), while adding the second pair Eqs. (B3- 
B4) gives the conservation law at the (l) — (6) interface, 
— j^} = j^} — j ^). From these two equations we get 
jty — j^} = jQ , which is a conservation law across the 
entire barrier. Finally, by combining the Eqs. (B1-B3) 
above we get 


3 


(r) 


TlT2 7 -(0 
1-R 1 R 2 j ^’ 


(B5) 


from which we can read directly the overall transmission 
coefficient across the barrier, T to t = TiT 2 /(1 — R±R 2 ). 
For Ti = T 2 = T this yields 

Hot = 2 j’-i __ i > (B6) 


which was given in Eq. (32) of the main text. Similarly, 
from Eqs. (B2-B4) we obtain 




R\ + R 2 — 2RiR 2 
1 — R\R 2 


di) 


(B7) 


from which we get the overall reflection coefficient across 
the barrier. R to t = (Ri + R 2 — 2 RiR 2 )/(1 — R\R 2 ). 


Appendix C: Finite bias 

The model presented in the paper does not include 
the small voltage bias used in the experiment to create 
a nonzero electric current between the left and the right 
contacts. A small finite bias voltage Vj, can be easily 
included, as detailed below. We assume that this bias 
is equally distributed, 14/2 between the left metal and 
the graphene and — 14/2 between the right metal and 
the graphene sheet. We also assume that the currents 
are small enough so that the equilibrium structure of the 
energy levels in graphene is not considerably changed. 
The bias only provides a slightly higher Fermi level in the 
metal from which electrons are injected at the left side 
and holes are injected at the right side. Due to electron- 
hole symmetry, the total current and the total noise can 
be calculated by adding the contributions corresponding 
to electron and hole transport . 

The electron injected from the left into the graphene 
sheet will have a Fermi energy shift in the metal con¬ 
tact region A£4 m + eX4/2, and a Fermi energy shift 
AEpc + e!4/2 in the suspended region. Similarly, for 
the hole injected from the right side, we will have the 
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shifts AE’fm — eI4/2 and A£fg — eVJ,/2. The current at 
a bias voltage 14 + dVJ, can be calculated in this model by 
taking into account that this voltage will be distributed 
onto both junctions, dj = a^dVb/2 + a^dVb/2. The 
measured conductance for a bias voltage 14 is defined as 
a = dj/dVb, yielding 


er = 




(Cl) 


Similarly, the differential noise for a bias voltage 14 is 
s = 1 [«(+) + s (-)j , (C2) 

and the Fano factor becomes 

F=~. (C3) 

(7 

In these equations, cr (±) and s^ 1 are the same as a 
and s introduced in Sec. II, but are now calculated for 
electrons with energies A£4 m ± e!4/2, AEfq ± e!4/2 
shifted from the equilibrium values calculated in Eqs. 


(10) and (14), and correspondingly shifted momenta (see 
Eqs. (22) and ( 29)), k = |A.Efm ± e!4/2| hvp and 
4 m = I^fm ± e 14/21 hv f- To summarize, 




s (±) = 


4e 2 IU r m ' n {^FM i4g } 


7T h 


dk T*' ± ' ) 

J. tot 5 


4e 2 !F i4g } 


7T h 


dk,. 


1 -T 


(±) 


(C4) 


TlotiW 


where is defined according to Eq. (32), but at en¬ 
ergies AEfm ± eVb/2, AEfq ± eVb/2 and at shifted mo¬ 
menta *4 m, 4g ■ 

Note that due to the assumptions above, the impact 
parameter po and the slope a will not change, since these 
are calculated from the equilibrium values eVo = AEfm~ 
AEfq, and a = v^I^oI/CfcI- Using this modified model 
we can account for the effect of a small bias voltage. The 
simulation show that, as expected, that this produces 
a slight broadening of the ideal (zero-bias) conductance 
and Fano factor, and does not have a significant effect 
for the determination of the doping parameters. 
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